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ABSTRACT 

We report the first detection of gas-phase CO2 emission in the star-forming 
region Cepheus A East, obtained by spectral line mapping of the V2 bending 
mode at 14.98 /xm with the Infrared Spectrograph (IRS) instrument onboard the 
Spitzer Space Telescope. The gaseous CO2 emission covers a region about 35" x 
25" in extent, and results from radiative pumping by 15 fim continuum photons 
emanating predominantly from the HW2 protostellar region. The gaseous CO2 
exhibits a temperature distribution ranging from 50 K to 200 K. A correlation 
between the gas-phase CO2 distribution and that of H2 S(2), a tracer of shock 
activity, indicates that the CO2 molecules originate in a cool post-shock gas 
component associated with the outflow powered by HW2. The presence of CO2 
ice absorption features at 15.20 toward this region and the lack of correlation 
between the IR continuum emission and the CO2 gas emission distribution further 
suggest that the gaseous CO2 molecules are mainly sputtered off grain mantles 

- by the passage of slow non-dissociative shocks with velocities of 15-30 km s~^ 

- rather than sublimated through grain heating. 

Subject headings: ISM: Molecules — ISM: Clouds, molecular processes — Star- 
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1. Introduction 

Cepheus A is a well-known site of star formation that has been studied extensively. 
Previous results have shown that it contains a series of heavily embedded far-infrared and 
radio-continuum sources, one of which dominates the luminosity of the entire region (HW2 
with 2.5x10'' Lq; Hughes & Wouterloot 1984; Evans et al. 1981). Ground- (e.g. Goetz 
et al. 1998) and space-based observations (e.g. Hartigan et al. 1996; Wright et al. 1996; 
van den Ancker et al. 2000) revealed that the multipolar outflow existing in this region 
consists of complex structures of shock-excited atomic and molecular gas. H2 emissions 
and line emissions from ionized species were detected along high-velocity (HV) and extra- 
high-velocity (EHV) outflows, indicating the presence of both dissociative (J-type) and non- 
dissociative (C-type) shocks, likely resulting from successive episodes of activity (Narayanan 
& Walker 1996). The protostellar object HW2 was determined to be the dominant powering 
source of the EHV outflow, while the source of the HV jet is still under debate (Goetz et 
al. 1998; Hiriart et al. 2004). The cavities carved into the surrounding molecular gas by the 
quadrupolar outflow are clearly seen in the NH3 maps obtained by Torrelles et al. (1993). 

In this Letter, we focus on new observations of Cepheus A East obtained with the In- 
frared Spectrograph (IRS) onboard the Spitzer Space Telescope {SST). The benefit of the 
.S'.S'T'data lies in the much greater sensitivity and the significantly better spatial resolution of- 
fered by the IRS compared to the Infrared Space Observatory (ISO). Such advances allow for 
the detection and spectral mapping of variations in the molecular and atomic line intensities 
on scales of a few arcsec instead of the previous few tens of arcsec or few arcmin scales, thus 
enabling our understanding of shock physics and chemistry at much finer scales. Detections 
of H2 S(l) through S(7), C2H2, [Ne II], [Ne III], [S I], [S III] and [Fe II] emissions as weU as 
absorption from CO2 and H2O ices will be discussed in a future paper. The present work 
will focus on the first detection of gas-phase CO2 emission toward Cepheus A East. Section 
2 describes the observations and data analysis. Sections 3 & 4 discuss our new results as well 
as the spatial distributions of gas-phase CO2 and H2 S(2) in the context of shock chemistry 
and outfiow activity. 

2. Observations and data reduction 

Spectral maps of two overlapping 1' x I'-square fields were obtained toward Cepheus 
A East with the IRS instrument onboard the Spitzer Space Telescope (SST) as part of 
Guaranteed Time Observer (GTO) program 113. The Short- Low (SL, both orders), Short- 
High (SH) and Long- High (LH) modules allowed for wavelength coverage from 5.2 to 25 
fj,m. Most of the data longward of 25 //m suffer severe detector saturation as a result of 
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very strong continuum emission and were not used. Continuous spatial coverage in the 
overlapping fields was obtained by stepping the slit perpendicular and parrallel to its long 
axis in steps of one-half its width and 4/5 its length, respectively. 

The data were processed at the Spitzer Science Center (SSC) using version 12 of the 
IRS pipeline. The Spectroscopy Modeling Analysis and Reduction Tool (SMART) software 
(Higdon et al. 2004) was used to extract wavelength-calibrated spectra at each spatial 
position in the maps. We used a set of locally-developed IDL routines specifically designed 
to remove bad pixels in the high resolution modules, to calibrate the fluxes for extended 
sources, and to generate spectral line maps from the extracted spectra (Neufeld et al. 2006). 

3. Results 

Figure 1 shows two examples of spectral line intensity maps we obtained using the SH 
data alone to allow direct comparison of the H2 S(2) and gaseous CO2 spatial distributions. 
The right panel displays the spatial distribution of the H2 S(2) intensity (12.28 /xm), while 
the left panel shows the intensity distribution for gas-phase CO2 (Q-branch of the V2 bending 
mode at 14.98 /xm). The H2 and the CO2 emissions mainly arise at the surfaces of the NH3- 
free cavities carved by the outflows. These emissions thus seem to result from interactions 
between the EHV outflow and the quiescent molecular gas traced by NH3 (Gomez et al. 
1999). Figure 2 (left panel) displays examples of summed spectra in the wavelength range 
relevant for this study. 

To constrain the physical conditions in the C02-containing gas, we generated synthetic 
profiles of the V2 CO2 band for temperatures ranging from 50 K to 900 K and compared them 
with the observed spectra (Fig. 2, right panel). We find that temperatures between 50 and 
200 K best fit the observed CO2 Q-branch emission (14.98 //m) over the entire C02-emitting 
region. Our observations, therefore, reveal the presence of a post-shock gas component much 
cooler than that measured with H2 pure rotational line transitions (e.g. van den Anckcr et 
al. 2000, T ~ 730 K from H2 S(l)-S(5) over the ISO Short Wavelength Spectrometer beam 
sizes 14"x20" and 14"x27"). 

To determine the column density associated with the gas-phase CO2 emission, one needs 
to identify the excitation mechanism. Prior studies (e.g., Gonzalez-Alfonso & Cernicharo 
1999; Boonman et al. 2003) showed that gaseous CO2 molecules can be excited: 1) by 
collisions with H and H2; 2) by radiative pumping into the 4.27 jiva. band; 3) by radiative 
pumping due to 15 jim. continuum photons emitted by dust local either to the CO2 gas 
component or to HW2. 
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For collisions to dominate, densities in excess of n ~ 10^ cm~^ would be required. How- 
ever, prior observations indicate densities from few~10^ to few x 10'' cm^'^ in the outflow 
region (Codella et al. 2005). Radiative pumping into the 4.27 /im band and subsequent cas- 
cade can be disregarded because the relaxation to the 15 /im band would produce detectable 
emission at 13.9 //m and 16.2 ^m that we do not observe. Resonant scattering by 15 yum 
continuum photons is, therefore, the most likely mechanism to excite the observed gaseous 
CO2 molecules. Continuum maps obtained from our data also indicate that the radiation 
field is dominated by dust emission close to HW2 - from a reflection nebula (Casement & 
McLean 1996; Goetz et al. 1998, Martin-Pintado et al. 2005)- rather than by dust emission 
local to the CO2 gas component. We, hence, conclude that the gas-phase CO2 molecules 
detected toward Cepheus A East are predominantly excited via radiative pumping by 15 //m 
continuum photons emanating from the protostellar region HW2. Resonant scattering in the 
CO and H2O vibrational bands, through continuum photons emanating from lRc2/BN, was 
also found toward the shock region Orion Peak 1/2 (Gonzalez- Alfonso et al. 2002). 

To compute iV(C02), we adopted the total luminosity and source effective temperature 
from Lenzen et al. (1984), as well as a distance of 690 pc for the Cepheus A East region, 
and we assumed that the projected distance of each spatial position is equal to its true 
distance from HW2. Under that assumption, the inferred CO2 column density is proportional 
to the measured intensity and to the square of the angular separation from the exciting 
source HW2. We also corrected the CO2 intensity measurements for extinction caused by 
the CO2 ice absorption. A comparison of the N{C02) map (Fig. 3, left panel) with the 
intensity distribution of the H2 S(2) hne (/(H2) S(2), Fig. 1, right panel) indicates that 
- like H2 S(2) - CO2 peaks at the NE position, coincident with the NH3 bridge between 
Ccp-A2 and Ccp-A3 (Torrelles et al. 1993) and at the HW5/HW6 positions. These striking 
similarities again strongly suggest that the CO2 emission arises in post-shock gas, traced 
by the warm H2 emission, and results from interactions between the EHV outflow and the 
ambient medium. The right panel of Fig. 3 displays the derived, extinction corrected, column 
density measurements of gas-phase CO2 versus H2 S(2) for spatial positions where the CO2 
line intensity (Q-branch) is greater than 0.75x10"^ erg s~^ cm~^ sr~^. Our estimate of 
the Spearman coefficient (p ~ 0.42 and significance level <1%) indicates that the column 
densities of the two species tend to increase together with spatial position. 



4. Discussion 



The presence of CO2 ice absorption features at most spatial positions in the gaseous 
C02-emitting region, their progressive weakening along the EHV outflow axis, as well as 
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the localization of the gaseous CO2 emission into a cool post-shock component, all pose the 
question of the origin of the gas-phase CO2 molecules in Cepheus A East. Was CO2 injected 
into the gas-phase by sublimation from C02-rich dust mantles after substantial grain heating 
or by grain sputtering caused by the passage of slow non-dissociative shocks known to exist 
in this active region? 

We searched for evidence of a correlation between the dust continuum emission, a tracer 
of grain heating, and the gaseous CO2 distribution in the C02-emitting region by comparing 
the CO2 column density measurements with measurements of the dust continuum emission 
at 7.3 fim, 14.8 yum and 18.0 /im. Wc found no obvious relation between A^(C02) and 
the dust continuum emission at these wavelengths. Furthermore, we determined that the 
dust temperature is never greater than the CO2 sublimation temperature (~ 90 K) at the 
locations of the CO2 column peaks. While grain heating likely occurs in the active region 
close to HW2, C02-rich grain mantle subhmation is not the dominant source of gaseous CO2 
in Cepheus A East. Sputtering off C02-rich ice mantles is, therefore, the most likely scenario 
in the present case. Our measurements of the CO2 ice column density range from ~ 10^^ to 
~ few X 10^^ cm~^ over the region exhibiting CO2 gas emission, indicating that only a few 
percent of the material has passed through mantle-destroying shocks. 

To derive the abundance of gaseous CO2, an estimate of A^(H2) in the CO2 gas com- 
ponent is needed. Because the gas-phase CO2 molecules arise from a post-shock component 
much cooler than that containing the H2 gas we detected, we cannot directly measure the 
total column density of the gas at T = 50-200 K from our data. Our abundance estimate 
will, hence, rely on previous indirect measures. Gomez et al. (1999) derived iV(H2) = 1.5 
X 10^^ cm~^ in the outflow using HCO+ emission measurements; a similar estimate was 
obtained by van den Ancker (2000) using far-infrared CO emission at the ISO/Long Wave- 
length Spectrometer angular resolution (75" beam size). Adopting this H2 column density 
yields an average gas-phase CO2 abundance of few xlO~^. 

This value should be regarded as a lower limit since the A'^(C02) might be underesti- 
mated if the true distance from HW2 is greater than the projected distance we assumed. A 
small fraction of the gaseous CO2 might arise from CO2 ice subhmation since some CO2 ice 
profiles, especially close to HW2, show substructures characteristic of grain heating (Ger- 
akines et al. 1999). A fraction of the gas-phase CO2 emission might also come from quiescent 
gas associated with Cep-A2 (knot at [R.A., Dec]= [7", 18"] in Fig. 1). While these additional 
sources of gaseous CO2 certainly contribute to the scatter seen in the right panel of Fig. 3, 
none of these sources can account for the correlation we observe between iV(C02) and A^(H2) 
S(2) over such a large spatial extent. 

The low gas-phase CO2 abundances detected by ISO (~10~^ instead of the expected 
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fewxlO"^) toward active star- forming regions with large CO2 ices reservoirs (e.g. van 
Dishoeck et al. 1996) fostered theoretical efforts to investigate the effects that shocks might 
have on CO2 chemistry. Charnley & Kaufman (2000) showed that for shocks faster than a 
critical shock speed, which is a function of preshock gas density, CO2 molecules sputtered 
off grain mantles are expected to be efficiently destroyed in the shock through reactions with 
H and H2. On the other hand, a shock speed higher than 15 km s~^ is necessary to allow 
efficient sputtering. 

Comparison of our results with the Charnley & Kaufman (2000) shock-model predictions 
indicates that slow C-type shocks with speeds of at least ~ 15 km s"^ interacted with the 
ambient gas over the extent of the C02-emitting region (see Fig 1). Adopting a pre-shock 
density of uh ~ 10^ cm~^ (Goetz et al. 1998), our abundance limit further indicates that 
the shock speed is no greater than 30 km s~^ since, at this density, higher shock speeds 
would efficiently remove CO2 from the gas phase, mainly through reactions with hydrogen 
atoms (see Fig. 3 of Charnley & Kaufman 2000). The CO2 emission might arise at the 
flanks of the bow shock modeled by Froebrich ct al. (2002). The average gas phase CO2 
abundance derived above (few xlO^'') is consistent with a scenario in which a few percent of 
the grain material has been subject to shocks fast enough to sputter CO2 ice mantles (with 
C02(ice)/H2 - few xlO'^). 
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Fig. 1. — Intensity distribution of gas-phase CO2 emission (Q-branch at 14.98 /xm) and H2 
S(2) pure rotational line emission (at 12.28 /xm) toward Cepheus A East is shown as black 
contours. The lowest intensity contours are 0.75 and 1.0 for CO2 and 1.0 for H2 S(2) with 
steps of 0.5 X 10~^ erg cm~^ sr~^. The (x) indicate the position of the radio-continuum 
sources HW2, HW4, HW5, HW6 and the EHV outflow edge (NE: AacosS = 27"; A5 = 
22"). The coordinates are offsets in R.A. (AacosS) and declination {A5) in arcscc with 
respect to HW2 (J2000: a =22h56ml7s.9 and 6 = +62°01'49"; Hughes & Wouterloot 1984). 
The grayscale shows the distribution of NH3 (1,1), a tracer of the cold quiescent molecular 
gas. The lowest contours arc 10 and 25 with steps of 25 mJy km s~^ beam~^ (Torrefies et 
al. 1993). The black dashed line delineates the mapped region in Short-High. 
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Fig. 2. — left panel: Summed IRS spectra of gas-phase CO2, H2, C2H2, [Ne II] and [Ne III] 
emission and CO2 ice absorption toward Cepheus A East. Individual spectra were summed 
over regions ~6"x8" in size, centered on the radio-continuum sources indicated in Fig. 1. 
Note the appearance of the gas-phase CO2 emission feature and the progressive decrease in 
the CO2 ice absorption feature when moving along the axis of the EHV outflow traced by the 
H2 emission and the NH3 cavities and away from HW2. Note also the presence of gaseous 
CO2 emission alone at the NE position, the farthest from HW2 in our maps. Right panel: 
Continuum-subtracted summed spectrum from the HW5/HW6 position. Overimposed (red 
trace) is the calculated CO2 band spectrum for T = 50, 200 and 900K. Note the progressive 
shift of the bandhead and the Q-branch widening with increasing temperature. 
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Fig. 3. — Left panel: Distribution of gas-phase CO2 column density, 7V(C02). The lowest 
contours are 1.5 and 3 with steps of 0.5x10^^ cm~^. The symbols and designations are as in 
Fig. 1. Right panel: N{C02) vs A^(H2) J = 4 (S(2) line) for the spatial positions exhibiting 
gaseous CO2 (Q-branch) intensities starting at 0.75 xlO~^ erg s~^ cm~^ sr~^. 



